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JiLiI Abstract 

We study to what extent the three-dimensional SU(A^) + adjoint Higgs theory can be 
1-^ ' used as an effective theory for finite temperature SU(A^) gauge theory, with = 2, 3. 

The parameters of the 3d theory are computed in 2-loop perturbation theory in terms 
^ . of T/Aj^, N, Nf. The perturbative effective potential of the 3d theory is computed to 

two loops for N = 2. While the Z(A^) symmetry probably driving the 4d confinement- 
deconfinement phase transition (for Nf = 0) is not explicit in the effective Lagrangian, 
it is partly reinstated by radiative effects in the 3d theory. Lattice simulations in the 
3d theory are carried out for N = 2, and the static screening masses relevant for the 
high-temperature phase of the 4d theory are measured. In particular, we measure 
non-perturbatively the 0{g^T) correction to the Debye screening mass. We find that 
non-perturbative effects are much larger in the SU(2) + adjoint Higgs theory than in 
the SU(2) -|- fundamental Higgs theory. 
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1 Introduction 



Corresponding to the known forces of nature there are two important finite temper- 
ature phase transitions in elementary particle matter: the QCD phase transition at 
T = O(IOO) MeV and the EW (electroweak) phase transition at T = O(IOO) GeV. The 
former has been studied intensely with numerical lattice Monte Carlo simulations 
but, due to the difficulties associated with treating dynamical quarks, conclusive state- 
ments about the properties of the physical QCD transition cannot yet be made. In 
contrast, for the EW case the problem can be essentially solved [§] by a combination 
of analytic and numerical means: first deriving by a perturbative computation 0|-|T^ 
a 3d effective theory Scs for the full 4d finite T theory ( "dimensional reduction" ) and 
then solving this confining non-perturbative 3d theory by numerical means |rT|-|TB|. 
The effective field theory approach has been intensively used for computations in high 
temperature QCD |T3-|T^, as well. 

In fact, there has been no doubt that dimensional reduction of QCD would work 
well at very high temperatures in the QCD plasma phase, T ^ T^, in the sense of the 
3d theory giving correctly the static correlation functions of the theory. The situation 
is quite different in the phase transition region J. Tc'. the effective finite temperature 
gauge coupling g'^{T) is becoming large so that the mass hierarchy g'^T,gT -C vrT 
needed for the construction of a simple local effective theory is lost. In particular, 
quarks play dynamically a crucial role and it seems that an effective theory using 
constant (in imaginary time) field configurations as essential degrees of freedom cannot 
be accurate (quarks are antiperiodic and are thus integrated out in this effective theory, 
their effect appearing only in the 3d parameters). However, it still seems well motivated 
to study how far one can actually go with dimensionally reduced QCD towards the 
phase transition region and the purpose of this paper is to do this. In comparison with 



earlier work [20|-|23| we go further firstly by determining the two continuum parameters 
of the effective theory in terms of T/Aqcd, N, Nf (quark masses are neglected) using 2- 
loop perturbation theory and the techniques developed in 0. Secondly, we are now able 
to extrapolate the lattice results to the continuum limit, using the lattice-continuum 
relations derived in 



n], H (see also |5 



Our conclusions will thus be different from 
those of the previous lattice studies [^-[|22[ for SU(2). 

Quite independent of its use as a finite T effective theory, the 3d SU(2)+adjoint 
Higgs theory is interesting because it has monopoles 00, E7[ which remove the photon 



from the physical spectrum and replace it by a pseudoscalar with the mass ||28[ 
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where mw is the perturbative mass of the vector excitation in the broken phase (in a 
gauge invariant analysis does not correspond to a physical state). A measurement 
of the photon mass with gauge invariant operators will permit one to make statements 
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about the monopole density in the broken and symmetric phases. These questions have 
recently been addressed in p9| . 

The paper is organized as follows. In Sec. 2 we discuss the derivation of the effective 
3d theory. In Sec. 3 we perform some perturbative estimates within that theory. In 
Sec. 4 we address the role of the Z(iV) symmetry in the effective 3d theory, and in 
Sec. 5 we consider the 3d lattice results. The implications of the analysis for the 4d 
finite temperature gauge theory are in Sec. 6, and the conclusions in Sec. 7. 



2 Relating the 4d and 3d theories 

Finite T QCD (with colours for the moment) is defined by the action 

S = f^dr j d'x^^F^^P;^ + Y.ni,D^{A) + m#,}, (2.1) 

where 

F^u = d,K - d.A^, - gfabcAlAl, (2.2) 
D^ = d^ + zgA„ A^ = AIT% (2.3) 

and where the gluon field is periodic (quark field antiperiodic) in r with period [3 = 1/T. 
In the matrix representation 

F^u = d^A, - d,A^ + %g\A^, A,] = {ig)-^[D^, D,]. (2.4) 

After renormalisation (we use the MS scheme), g"^ becomes scale dependent, and to 
1-loop 

Some useful group theoretical relations for the SU(N) generators are given in Ap- 
pendix A. 

Upon dimensional reduction, the field Aq becomes an adjoint Higgs field and the 
general form of the super-renormalizable Lagrangian of the 3d effective theory can be 
written down: 

Leff[4^Ag] = ^Ft^Ft^+TT[D,,Ao][D„Ao] 

+ mlTTAl + XA{TiAiy + XA[TTAl-^{TTAiy]. (2.6) 

The reduction process does not generate terms proportional to TrAg. Operators of 
higher dimensions are parametrically of the form g^A^ and can be neglected relative to 
the retained term g'^AQ as long as Aq <^ 2tiT / g (their contributions to the correlators 
we will measure are also of higher order) . 
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Note that at very high temperatures, it is even possible to integrate out the Aq- 
field [|, |, [1^ . Here we want to go as close to Tc as possible, and hence we keep Aq 



in the effective Lagrangian. 

For general N two independent quartic couplings appear in eq. ( ^.6[ ). We shall take 
in this paper = 2, 3 for which eq. ( [A.4|) implies that one can take = 0. The case 



= 5 with both couplings is treated in |^ 

The effective theory thus depends on three dimensionful couplings: g^, m|,, A^- 
Instead of these one can use one dimensionful scale and two dimensionless parameters, 
chosen as 

al y^'^, (2.7) 

The 3d theory is super- renormalizable and g^, \a are renormalisation scale independent 
while m|, is of the form 

ho = 2{N^ + l){Ngl-\A)\A, (2.8) 



where is a constant specifying the theory. The scale in the definition of y in eq. 
is chosen as the natural one, g^. It is noteworthy that there is no (/g-term in ( |2.8|) . 

The process of dimensional reduction now implies finding the relation between the 
physical parameters of finite temperature QCD and gl,y,x. The Lagrangian parame- 
ters of QCD are after renormalisation Aj^, mj(/i); the physical parameters are hadron 
masses (one mass to set the scale and the rest as dimensionless mass ratios) . As hadron 
masses are entirely non-perturbative, it is conventional to use as(m^) and define Aj^ 
as the scale for which as diverges when evolved to smaller scales. 

We will derive the parameters of the 3d theory so that the relative errors are of the 
order 0{g'^). This requires a 1-loop calculation for the gauge coupling, but for the mass 
parameter m|, and the scalar coupling A^ one needs a 2-loop derivation. The actual 
reliability of this calculation is to be discussed below. 

To 1-loop (tree level for g^) the calculation gives ||, |^ 

gl = g\f^)T, (2.9) 

\{N + lNf)g'{f^)T\ (2.10) 



.g\f^)T 



Xa = + (2.11) 

Xa = (N-Nff-^. (2.12) 

We take Nf flavours of massless quarks, although the dependence on mass thresholds 
could also, in principle, be included. 
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At this level the scale is unspecified; thus a 2-loop derivation of the effective theory, 
constituting a resummation of the perturbative series and establishing the scale at 
which the 3d parameters are to be evaluated, is needed. Indeed, the 3d couplings are 
scale independent (note that the 3d scale dependence in eq. ( ^.81 ) is of order and 
thus does not yet enter at this level of the 4d— i>3d reduction) so that the /i dependence 
of the 2-loop result must be the following (we only discuss = 2, 3 so that Aa is 
irrelevant): 



9l 
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where 



L 



22N ^ a 
—-In — 

3 ht 



ANf , 4u 
3 ht 



(2.13) 
(2.14) 
(2.15) 

(2.16) 



and the Cj are constants to be found. 

The derivation of the parameters can be most easily made using the background field 
gauge. Calculating to 1-loop the contribution from the n 7^ modes to the zero mode 
correlators, one gets the relations between the 3d and 4d fields: 



where ^ is the gauge parameter and 



3 V /iT J 3 \ fiT 

.^f221„A + l) + ^f41„l^i 

3 V fix J 3 V i^T 



(2.17) 



f2.18) 



is the standard thermal scale arising in perturbative reduction in the MS scheme 



The 3d gauge coupling can be read directly from the gauge independent A^A^ 

For the other parameters, one needs the 2-loop 



normalization factor in this gauge |17 



correlators at zero momenta of the Ag-fields in the 4d and 3d theories, to separate the 
contributions coming from the n 7^ modes. These can be most easily derived with the 
effective potential. The full 4d 1-loop effective potential can be read from eqs. (4.2) 

0: 



(with 3/4 4/3), (5.2) of ||3T[, and the 2-loop effective potential from eqs. (4.3), 
(5.4) of The 1-loop potential is gauge independent, but the 2-loop potential in 
these equations corresponds to ^ = 3, since the extra rescalings added in eqs. (3.17), 
(3.18) of 1^ vanish for that gauge parameter. For the mass parameter ml one has to 



subtract the contribution from 3d, which has to be calculated separately, but for Xa 
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there is no 3d contribution in the 4d 2-loop effective potential and thus the coefficient 
of the quartic term gives directly the n 7^ contribution. Finally, to get the 2-loop 
results for the 3d parameters, one needs the terms arising when the rescalings in (|2.17|) 
combine with the 1-loop results. Alternatively, the 2-loop mass parameter can be 
read directly from |19[. Note that the 2-loop cubic terms in the 4d effective potential 
are not used in the derivation of the 3d parameters and are not reproduced by the 3d 
theory (they are higher order contributions when SA<ttT). 
The computation described above gives 

N 

c, = y, (2.19) 



lON^ + 2N] + 9Nf/N 



(2.20) 



(2) 7/3 - 109iV./96 2 ,^ , , 

cf = V2 - 23A^,/18 2 

' 1 - Nf/9 3 ^ ^ ^ 

Note that for the computation of the free energy of QCD in the symmetric phase to 
order g^T"^ through a 3d theory [|T^, one only needs c^; the constants Cg and c[^'' give 
higher order contributions. 

Another useful representation of eqs. ( P^.13 )-( ^.15 ) can be derived by choosing dif- 



ferent renormalization scales for different parameters in a way that 1-loop corrections 
to the 3d parameters vanish. In general, the solution of L -|- c = is 

/-3c + 4A^f ln4\ 
^^ = ^^T exp 22iV-4Ar, J " ^2.23) 

Using dU) and (|1|) this gives (for Nj = 0): 

9l 247r2 



T llA^ln(6.742T/A^^)' 



(2.24) 



Secondly, for x one obtains 



gl llN-2Nf ln(/iT/i/AMs) ' ^ ' ^ 

where fi (= e^^/^^ for Nf = 0) is obtained from (|2.23|) with c = 2c[^^ —Cg. This equation 
gives quantitatively the relation between the 4d theory variables T/Ajj^, N, Nf and the 
3d effective theory variable x. For Nf = 0, = 2, 3 one has 

_ 6 + A^ 1 

~ TiaT ln(5.371T/Aj^) ' ^^'^^^ 
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Finally, concerning y = fnj){g1)/gt note that the RHS of eqs. (|2.13|) - (p.l5|) only 
depends on and Nf. Eliminating g^i^n) one obtains a //-independent relation 

between the dimensionless scale independent variables x. To leading order y = 
mjj/g^ ~ l/g'^;x= \a/ g1 ~ S'^ so that ?/ ~ 1/x. The complete relation is 

(^^3). . _ (9-iV,)(6 + iV,) 486 -33iV,-lliV|-2iV3 
For Nf = these have the simple forms 

- i + 1^5-8^(1 4^ + (2.30) 



With the use of known lattice data for the phase transition in pure gauge theories in 
4d for = 2, 3 one can define the value of x corresponding to the critical temperature. 
We have: TJA^ = 1.23(11) for SU(2) and TJA^ = 1.03(19) for SU(3) according 
to [^. Thus, for = 2 the value of x corresponding to Tc is about Xc = 0.2 {xc — 0.17 
for SU(2) and Xc — 0.14 for SU(3)). To have a feeling of the accuracy required for the 
assessment of the lattice results below, the leading value of y{0.2) is 0.1126 to which a 
22.5% 2-loop correction 0.0253 is to be added. One might then estimate that the next 
(omitted) term is roughly 5% so that the theoretical result is y{x = 0.2) = 0.138(6). 

Quite surprisingly, the power series defining the mapping of the 4d parameters to 
the 3d ones seems thus to be quite convergent even at the critical temperature. How- 
ever, this does not prove that the 3d theory can adequately describe the confinement- 
deconfinement phase transition at high temperatures. There is another important 
criterion for the applicability of dimensional reduction, namely that the typical 3d 
mass scale must be much smaller than vrT, since only then the integration out of the 
non-zero Matsubara modes is self-consistent. As we will see, it is this point which does 
not allow an effective 3d description near the critical point. 

Summarizing, to the extent that finite T QCD can be regarded as an SU(N) gauge 
theory with Nf massless quarks characterized by the physical quantities T/A^, N, Nf, 
a 3d effective theory given by the super-renormalizable Lagrangian ( p. 61 ) with the cou- 
plings g\,y,x (eqs. (p^) ) can be derived. The relation between these two sets is in 
eqs. ( p^ , ( piel) , ( p:27D and <^J2^ ). 



Note the crucial difference in comparison with SU(N) + fundamental Higgs theories, 
in which there is a Higgs potential with two parameters already at the tree level. Then 
(for N = 2) the variable x (the Higgs self-coupling in the effective theory) is essentially 
the zero temperature Higgs mass and y (the scalar mass in the effective theory) is 
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~ (T — Tc)/Tc. The whole x > plane corresponds to some physical SU(2)+Higgs 
finite T theories. For pure SU(2), in contrast, both x and y are given by T/Aj^, A^, Nf 
and for fixed iV, Nf only one curve yM-*?,d{,x) corresponds to a physical 4d theory. 
Assume the effective theory has a phase transition along some curve y = ydx), which 
we shall soon determine with lattice Monte Carlo simulations. Then this transition 
could correspond to a physical 4d transition only if it intersects the curve y4d-^3d{x) in 
a region where the derivation of the effective theory is reliable. 



3 Perturbation theory in 3d for N = 2 

The aim now is to study the phase structure of the 3d theory defined by the Lagrangian 
on the plane of its variables x, y. On the full quantum level the study has to be 
based on gauge invariant operators and will be carried out by numerical means in 
Sec. 5; only this gives reliable answers. However, it is also quite useful to fix the gauge 
and study the problem perturbatively. On the tree level the answer then is obvious: 
there is a symmetric phase for ?/ > (m|) > 0) at all x > and a broken phase for 
?/ < 0. To get a more accurate result, one shifts the field by 



Aa 
^0 



obtains 3=1+2 scalars with masses 



D 



+ 3A, 



^3(1/ + X0^), 



two massive vectors with mass 



(3.1) 



(3.2) 



2 J,2 



and one massless vector. The 1-loop potential in the Landau gauge is 

1 



loop 



127r 
1 



where 



127r 



1 



4my + m\ + 2m\ 
403 + + 3x02)3/2 ^ 2(y + a;02)3/2 



y(^) = y+ -(20x - lOx^) In 
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(3.3) 



(3.4) 



(3.5) 



IGtt^ /i 

Let us consider first the high temperature limit of the 1-loop effective potential. It 
corresponds to the case x — >■ 0. The leading terms are: 



Vi-\oop/gl 



1 



X02 



2 \2 



37rx 



+ 



X 



V dn'^xy 



(3.6) 
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Two degenerate states are obtained when the last term vanishes. Prom this one finds 
that 



0symm 0, ^broken 

23/2 

0-3 = 



8l7r3x5/3' 



(3.7) 
(3.8) 
(3.9) 



where is the 3d interface tension divided by g^. Further, the upper and lower 
metastability points y±{x) are given by 



X 



(3.10) 



The contributions to the 2-loop potential from vectors {-^), ghosts {^=) and scalars 
— ) are (without a factor l/(167r^)) 



89, 



= gln4[AH{mT, mr, 0) + H(mT, 0, 0) - — ] , 




CO = 



5f3mr-H"(mr, 0,0), 



2 2 



^^^^^ 




-293[^ssv{m2, m2, 0) + 2Dssv{mi, ma, mr)], 
2glmT{mi + m2), 

-g^(j)^[Dvvs{0, rriT, ma) + 2Dvvs{mT, mr, mi)], 
-A^0^[3i7(mi, mi, mi) + 2i7(mi, ma, ma)], 
-Ayi[3m^ + 4mima + Sma], 



(3.11) 
(3.12) 
(3.13) 
(3.14) 
(3.15) 
(3.16) 

(3.17) 
(3.18) 



where 



Dssv{mi,m2,M) = 

(M^ - 2mi - 2ml)H{M,mi,m2) + -^(mi + ma)[M2 + (mi - ma)^] 



mima 



+ 



1 



M2 



m^ - ma) [i/(M, mi, ma) - iJ(0, mi, ma)]. 



(3.19) 



Dvvs{Mi,M2,m3) = 

\h{M^, M2, mg) + — ^[(Mi - M^f - m^iM, + M2 + m^)] 

Z 4:IVllIVl2 



+J]^p]^{'^t[H{M,, M2, ma) - H{M,, m,, 0) - HiM^, m,, 0) + H{m3, 0, 0)] 
+ (Mf - 2M2m^)[i/(Mi, M2, ms) - H{Mi, m^, 0)] 

+ (M24-2M>2)[^(Mi,M2,m3) -/f(M2,m3,0)]|. (3.20) 



The sunset function is 



1 1 ' 

H {nil, 1712, ms) = —^[\n . . + 77)- (3-2i; 

loTT"' V mi + m2 + ma 2/ 



In (|3.11|) - (13.201) the factor l/(167r^) appearing in H should be suppressed. 



The use of 2-loop effective potential allows the construction of the critical curve ydx) 
beyond 1-loop level. The numerical computation is shown in Fig. |l|. 



4 Z(A^) symmetry in the 3d theory 



For Nf = 0, hot SU(A^) gauge theory exhibits a Z(A^) symmetry ^ so that there 
are equivalent ground states. In one of them Aq = 0, in the others Aq ~ 2ttT / g. 
In the weak coupling limit the field in the other minima thus becomes large. One can 
associate the QCD phase transition with the breaking of this symmetry so that at high 
temperatures one sits in one of the equivalent minima, say Aq = 0. As one goes to 
lower temperatures, the barrier between the minima is becoming lower and at T^, the 
symmetry is supposed to get restored. If one wants to apply a 3d theory down to 
temperatures close to Tc, it is thus important to discuss the role of the Z(A^) symmetry 
in 3d. 

The construction of the super-renormalizable effective 3d field theory in Sec. 2 re- 
quires small amplitudes of the adjoint field Aq, Aq ^ 2'kT/ g. In practice, this means 
that one works around the = minimum in the reduction step, so that the ef- 
fective theory describes reliably only fluctuations of magnitude 5 A < ttT around that 
minimum. Thus, the price one pays for the simple effective theory is that the Z(A^) 
symmetry is not reproduced by it. 

Some remnant of the Z(iV) symmetry nonetheless remains in the 3d theory. The 
effective theory will on its critical curve still contain metastable states but these 
are not completely equivalent, e.g., the correlators measured in them are different. As 
discussed, only the correlators in the symmetric phase {Aq = 0) reliably represent 4d 
physics. 
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Figure 1: The critical curve y = ydx) (multiplied by x) computed from 1- and 2-loop 
potentials in the SU(2)+adjoint Higgs theory. Two of the curves are plotted for the 
scale choice fi = g'^; fi for the second 2-loop curve is determined by optimization 0. 
The scale dependence should give an estimate of higher order corrections. The thin 
curve is yM^sdix) from eq. ( |2.29| ). For a; ^ all curves approach 2/(97r^) = 0.0225. 
The non-perturbative critical curve computed numerically is given in Fig. |^. 



Indeed, within the interval < < 27iT/g, eq. ( |3.6| ) at y^ is seen to be the same 
as the 4d 1-loop potential in an background: 



V{A 



rp4 



54(0) 
1 



2B. 



\27iTj 



90 3^ ^ °^ 



(4.1) 
27rT 



where 



BAx) 



30 



+ (xmodl)^[l — (xmodl)]^. 
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(4.2) 



The effective potential (|4.1| ) exhibits the Z(2) symmetry of the full 4d theory as the 
periodicity in A^ with period 2'kT / g\ this symmetry is not explicit in the effective theory 
on the tree level, but it is partly reinstated by radiative effects as seen in eq. (p.6|). This 
is not surprising: in deriving the 4d 1-loop potential one performs a frequency sum over 
n = 0, ±1, ±2, .... In deriving the couplings of the 3d theory one performs basically a 
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frequency sum over n = ±1, ±2, and the n = mode enters when computing the 
effective potential within the 3d theory. Note also that ( |3.7| ) is exactly the same as the 
leading term in ( p.29| ). 

Thus, while the Z(2) symmetry is not explicit in the effective action, the second 
degenerate minimum is generated radiatively. Even with all higher order corrections in 
the non-perturbative solution of the 3d theory there will be two (for N = 2) degenerate 
minima at some values of the 3d parameters, but this will no longer take place at any 
temperature T > Tc. Moreover, while these states are completely equivalent in the 
original 4d theory, the correlators measured in them will be different in the effective 
theory, the symmetric phase being the physical one. For SU(3) the full 4d theory has 
3 equivalent states while in the 3d theory presumably only two of them have the same 
correlators; these are the unphysical ones corresponding to a large value ~ 2t:T / g of 
the fields. 



5 Lattice analysis 



5.1 Discretization 

The lattice action corresponding to the continuum theory ( p.6| ) is S" = Sw + Sa where 



'^^ = /3GEE(l-^TrP,, 

X i<j ^ 



(5.1) 



is the standard SU(2) Wilson action and, in continuum normalization (Aq = A!^a) 
and for Aa = 0, 



E2 



where 



HxaAlix) - Tr aAo(x)t/i(a;)v4o(x + i)U\[x) 
+ E (^DCf)^Tray4Q(x) + aAA(Tr aA, 

X 

4 



(5.2) 



(5.3) 



is the lattice spacing and is the bare mass in the lattice scheme. Renormalisation 
is carried out so that the physical results are the same as in the MS scheme with the 
renormalized mass parameter 



2 / 20x-10x2 gl\ 4 

m.(^) = (y + ^^ln-j^3. 



The resulting counterterms have been computed in |23] with the result 



(2O2; - 10x2)(ln — + 0.09) +8.7+ 11.6s 



(5.4) 



(5.5) 
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where S = 3.1759114. There are no higher order corrections to this relation in the 
hmit a — i> 0. For 1-loop 0(a)-improvement, see p5[ . 

The action can be written in different forms by rescahngs of Aq. For example, 



replacing Aq by an anti-Hermitian matrix 
of the form 



SA = f3Aj2 ^Tr Ao{x)U,{x)Ao{x + i)Uj{x) + J2 



by aAl = — /S^^o/^' gives an action 

(5.6) 



-P2lTrAl + /3,i^TTAir 



where 

/32 _ 8y S(4 + 5x) 
Pa Ph It^Pg 



(20a; - 10x2)(ln -pc + 0.09) + 8.7 + 11. 6x 



and 



X 

Jg 



(5.7) 



(5.8) 



Since there are two dimensionless parameters, one more arbitrary choice is possible. 
In fundamental Higgs theories one often scales the coefficient of the quadratic term to 
unity; here we choose Pa = Pg- 

As the final relation, one wants to determine the continuum theory value of (AgAg) 
when doing simulations with the lattice action (|5?T|) + (|5.6| ) with Pa = Pg- The answer 
is 

{A-A-) P'g^ 3S 3 



^2 °^ AtiPg 



Ph' 



log ^ + 0.67 



(5.9) 



5.2 Simulations 

The aim of the simulations is as follows: 

• Find the critical curve y = ydx) and, in particular, its endpoint Xend. This is 
done by measuring the distributions of Tr A^ (or some other operator) at fixed x, 
finding a two-state signal and the V ^ oo,a limit of the value of y at which 
this happens, 

• Measure the correlator masses on both sides of the transition line and in the 
cross-over region x > Xend- 

It is important to estimate the required lattice volume V and constant a (or Pg)- In 
the broken phase, on the tree level we have at least two relevant mass scales, the adjoint 
scalar mass ~ gT from eq. (|3.2|) and the vector mass (7(^0) = 27rT from eq. ( p.3|) . Thus, 
in order to describe accurately the broken phase, we must demand 



a«^«^«Na. 



(5.10) 
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These convert to 

/?G»;r^, iV»^/5Gv/S. (5.11) 
Sirx 4 

For the symmetric phase the scale 27tT is absent, so that the requirement on (3g is not 
quite as demanding, /^g ^ 1. 

The simulations were performed with a Cray C90 at the Finnish Center for Scientific 
Computing, and the total cpu-time usage was 760 cpu-hours. 



5.3 The phase diagram 

We locate the transition line y = ydx) by using Tr^Q as an order parameter. For each 
lattice size, we define the pseudocritical y with the following commonly used methods 
3^ : (a) maximum of the susceptibility xi^o) = ((TrAg — (TrAg))^); (b) minimum 



of the 4th order Binder cumulant Bl{AI) = 1 - {{Tr Al)'^) /{3{{Tr ADY); and, when 
the transition is strongly 1st order, (c) the "equal weight" criterion for the probability 
distribution p(Tr Aq). For any finite volume, these criteria give different pseudocritical 
values of y, but they all extrapolate to the same limit when V — > oo. This extrapolation 
is shown in Fig. ^ for x = 0.20. 

The analysis is repeated for several corresponding to different lattice spacings 
through eq. (p.3|). In our case, as long as the conditions ( p.llj ) are satisfied, we found 



no appreciable lattice spacing dependence in the results. This is also evident in Fig.|^. 

Our result for the phase diagram of the continuum SU(2)+adjoint Higgs model is 
shown in Fig. |^. As already pointed out in the phase diagram consists of a first- 
order line which terminates, so that the two "phases" are analytically connected. We 
find the endpoint to be close to Xend = 0.3. As can be seen from the figure, the results 
are independent of /3g well within the statistical accuracy. However, the infinite volume 
extrapolation is essential, as shown by some large but finite volume points in Fig. |^. 

The thick transition line ydx) in Fig. |^ has been obtained with a fit to the oo 
-extrapolated results, as 

xyc = 2/(97r^)(l + 9/8x + Ax^^^ + Bx^ + Cx^^^ + Dx^), (5.12) 

where A = 45.1(2.9), B = -214.7(20.2), C = 350.7(44.5), D = -206.7(31.4). The fit 
has been constrained to join the curve xy^d^^dix) and its slope when x —>■ 0. We used 
fractional powers of x in the fit since in the perturbative expansion of the critical curve 
terms like g'^/mo ~ y/x appear. There can also be logarithms in the expansion. 

It is interesting that the behaviour of the 3d SU(2) + adjoint Higgs model is quite 
different from that of the 3d SU(2) + fundamental Higgs theory studied earlier |12| 



First, in the fundamental case the value of yc was smaller than the 2-loop perturbative 
result for all values of x, while in the adjoint case the situation is opposite for x < 0.2 
(compare Figs. |l|, |^). Second, the magnitude of the higher order perturbative or non- 
perturbative effects is much larger in the adjoint case. To get a feeling of the difference. 



13 



X = 0.20 




0.180 



0.170 



0.160 



0.150 



0.140 

0.000 0.005 0.010 0.015 0.020 



Figure 2: The Uc extrapolation to 1/ — oo for x = 0.20. Within the statistical errors, 
there is no lattice spacing Pg dependence. 

suppose that the deviation from the perturbative result comes entirely from a non- 
perturbative energy shift ai Aq = 



(5.13) 



where Ap is some constant to be determined numerically. For the fundamental case 
this was estimated to be ~ —0.08 at x ~ 0.06. A first order estimate of the change 
in the critical value, 

" -5V, (5.14) 



^SmlAl 



gives in the present case 



8 Sy, 



Apc^- — ^. (5.15) 

For example, for x = 0.04 we have Ap ~ 10 (~ 100 times larger than in the fundamental 
case!) and for x = 0.12, Ap ~ 0.4. We do not have any clear understanding of this huge 
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Figure 3: Data points in the limit ^ oo, a — for the critical curve y = ydx) 
(multiplied by x). The thick line is the 4th order fit ( 5.12|) to the = oo extrapolated 
data. The dashed line marks the region where the transition turns into a cross-over. 
The straight lines are the 4d— >-3d curves of eq. ( p.27| ) marked by the value of Nj. The 



top scale shows the values of T/Aj^ corresponding to the values of x for = 2, Nf = 0, 



obtained using eq. ( p.26| ). The physical implications of the figure are discussed in Sec. 6. 



difference. From the point of view of perturbation theory, both cases are quite similar. 
The only difference we see is that the adjoint theory contains instantons (monopoles) 
in the broken phase, while the fundamental does not. 

In Fig. ^ we show the probability distributions p{Aq), measured with various values 
of X along the transition line ydx). When x ~ 0.04, the transition is extremely strongly 
first order, as evidenced by the strong separation of the two peaks in the distribution. 
Indeed, for x < 0.12 we use the multicanonical algorithm in order to enhance the 
tunnelling from one phase to another - this is absolutely necessary in order to ensure 
the correct statistical weights for the phases. When x increases, the transition becomes 
weaker, and at x = 0.40 we see no sign of a transition any more. 

Due to the strong 1st order transition, there is a substantial metastability range 
in y around the transition line at small x. Within this range the thermodynamically 
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unstable state can be stable enough so that, for practical purposes (~ measurements), 
it behaves as if it were stable. In Fig. ^ we show the lower end of the symmetric phase 
metastability range for x = 0.04, 0.12 and 0.20. This has been obtained by reweighting 
the ydx) -distributions (Fig. ^) to smaller values of y until the peak of the distribution 
corresponding to the symmetric phase vanishes. 

The behaviour of (AqAq) /2g^ when crossing the transition line or its continuation is 
shown in Fig. ^. 



5.4 Correlator masses 

The list of gauge-invariant composite operators of lowest dimensionality is: 

• Dim = 1: the scalar Tr Aq, 

• Dim = 3/2: the vector hi = eijkAQF^^, 

• Dim = 2: the scalar (TrAg)^, 

• Dim = 3: the scalars (Tr A^)^ and G = and the tensor 



Perturbatively the Lagrangian ( p.6| ) for = 2 in the symmetric phase has 3x2 
massless vector (A'}) and 3 scalar (Aq) degrees of freedom. Because of confinement, 
only SU(2) singlets appear in the spectrum. Thus, the above mentioned operators 
represent bound states in the symmetric phase. For example, TrAg is a bound state 
of two adjoint scalar "quarks", hi is a bound state of an adjoint quark and light glue, 
and G is a scalar glueball state. 

In the limit of large 3> 1 the SU(2) theory under consideration in the symmetric 
phase can be considered an analog of QCD with heavy scalar adjoint quarks. Thus, a 
number of statements can be made on the expected spectrum with the use of the heavy 
quark expansion. First, the lowest energy states are different types of glueballs with 
niG ~ Qs- Next, we have a bound state hi of the scalar quark and gluon with the mass 
TUfi = rrii:) + 0{g1), mj:) = y/yg^. The leading correction to rrih can be determined: 

mh-mD = T^\og^/y + —gl + ... (5.16) 
ZTT Zn 

The coefficient in front of the log is a result of a perturbative computation of the pole 
Aq mass in whereas the constant term proportional to (?| cannot be computed 



perturbatively. In fact, the mass of this bound state coincides with the non-perturbative 
definition of the Debye screening mass, discussed in Finally, the mass of the bound 
state is expected to be 

m{Al) = 2^gl + 0{gl). (5.17) 
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Figure 4: The probability distributions of AQA'^/2gl = -^J^x ^o^of^Qi: measured at 
various locations on the transition line y = yc{x). 



In the broken phase the perturbative spectrum consists of 2x3 massive vectors iW^), 
mw — 25'|/(37ra:), 2 massless vectors (the photon 7) and one scalar degree of freedom 
(the Higgs) with mass ~ m^. The 2x3 massive vector states have U(l) charges 
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Figure 5: The behaviour of {A'^Aq) /2g'^ in continuum normahzation (see eq. ( |5.9| )) 
when crossing the phase transition at x = 0.12 and the cross-over at x = 0.35. 



and thus are not gauge invariant. Moreover, in 3d the U(l) Coulomb potential is loga- 
rithmically divergent, and thus the U(l) charged states (like W) must be absent from 
the spectrum. Therefore, the following set of gauge-invariant operators should describe 
the spectrum of the theory in the "broken" phase: Tr Aq for the Higgs excitation, hi 
for the "photon", and G for the bound state of the massive W^s. In fact, to obtain in 
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the unitary gauge Aq = A^/ AqAq = 5^3 a pure photon operator, one may take 

7,, = - -eatcA^,{D,Ao)\D,Aor. (5.18) 

fl'3 

However, in our measurements this operator was dominated by the first term, yielding 
similar results as the operator h^. As discussed in the introduction, the broken phase 
tree-level massless photon is actually replaced by a massive pseudoscalar excitation. 

In the following we constrain ourselves to the study of three different operators, 
Tr Aq, hi and G. Our strategy is the following: first, we are going to consider correlation 
functions near the 3d transition line. This will allow us to compare masses of different 
excitations in the vicinity of the phase transition, as well as to determine correlation 
lengths in the symmetric phase (as we will discuss later, this is the physical phase from 
the 4d point of view) along the 4d 3d dimensional reduction curve. Second, we make 
measurements at some fixed x for different values of y in the symmetric phase, in order 
to check different hypotheses concerning the spectrum of the bound states and in order 
to determine non-perturbatively the Debye screening mass and the 0{g^T) correction 
to it. 

The correlation functions are measured in the direction of the X3 -axis, and to enhance 
the projection to the ground states, we use blocking in the (xi, X2)-plane. The fields 
are recursively mapped from blocking level (k) — >■ {k + 1), so that the fields on the 
{k + l)-level lattice are defined only on even points of the (/c)-level lattice on the 
(xi, a;2)-planes: 



Ai'^'\y) = -Ai'\x) + ^ E ut\x)Ai'\x + t)ut^\x) (5.19) 

O LI O 



and {i = 1, 2,j 7^ i) 



Ur'\y) = f/f^(x)f/f)(x + .), (5.20) 

U:^'\x) = \ut\x) + \ E U^'\x)ui'\x + j)U^'^\x + t). (5.21) 

i=±i,2 

The scalar and vector operators above are then calculated from the blocked fields. 
The blocking is repeated up to 4 times, and the correlation functions are measured for 
each blocking level separately. For each operator and coupling constant we select the 
blocking level which yields the best exponential fits; typically level 3 or 4. 

Because we expect the vector operator hi to yield a very light mass value in the 
broken phase, we measure it in the lowest non-zero transverse momentum channel (for 
all of the blocking levels), in addition to the zero momentum channel: 

03(^3) = 4 E A^,F^,e^'-^^/''s _ (5.22) 
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Figure 6: The correlator masses of the operators Tr and hi = eijkTr AqFji. measured 
when crossing the phase transition at a; = 0.12 and the cross-over at a; = 0.35. 



The screening mass is then extracted from the asymptotic behaviour of {O3{0)Os{xs)) oc 
exp[— VTxs], where 

W = ^{2ti/NsY + m2. (5.23) 

The results for the masses both in the 1st order phase transition region (x = 0.12) 
and the cross-over region (x = 0.35) are shown in Fig. At the phase transition the 
masses exhibit a prominent discontinuity. In the broken phase y < tjc there is a photon 
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of very small mass and a scalar the mass of which is close to the mass parameter 



ttid; in these units y2/3gT/g'^T = y2/{9n'^x). In the broken phase only the non- 
zero transverse momentum operator gives statistically significant results for h^. The 
negative mass values in Fig. |^ notify negative m^-values, as obtained from the relation 
( p.23| ). The photon mass is statistically consistent with zero but, as discussed already 
in , a small non-zero photon mass is not excluded. This is required for the analytic 
connection between the two "phases" ; otherwise the photon mass could act as an order 
parameter separating the phases (this, in fact, is what happens in the U(l)+Higgs case, 
the Ginzburg-Landau theory |^). In the symmetric phase y > Uc there is a scalar 
Aq — Aq and a vector Ai — Aq bound state of rather large mass. In the cross-over region 
the mass of the scalar state dips a.t y = yc while that of the vector state increases 
monotonically. 

Now, consider the masses of A^ and hi in the symmetric phase along the dimensional 
reduction curve y = yM^3d{x), Fig. ^ These measurements are possible since for x < 
0.235 the symmetric phase is metastable with a sufficiently long tunneling time, while 
at X > 0.235 it is absolutely stable. As one can see the expectation m{Al)/m{h) ~ 2 
is strongly violated in the whole region of x studied, x > 0.05. This indicates that the 
corrections to the leading order results ( |5.16| ), ( |5.17| ) are large. In order to clarify this 



point, we study different correlators at large values of y and and fixed x, where the 
heavy quark expansion is expected to be valid. 

In Fig. P we show the dependence of the masses of the bound states on y at different 
values of X, X = 0.01,0.04,0.10. As expected, the glueball mass does practically not 
depend on y or x, while the Aq and hi masses have some dependence on x. Moreover, 
at large y, m{AQ) > m{h) as it should be. A fit to the data allows to fix the unknown 
coefficient c in ( p.l6| ). Only sufficiently large values of y should be used, ^/y > c, to 
make sure that the correction is smaller than the leading contribution. We found for 
y>2: 

X = 0.10 : c= 1.74(16) • (27r), 
X = 0.04 : c= 1.53(17) ■ (27r), 

X = 0.01 : c = 1.39(20) ■ (27r). (5.24) 

This indicates a slight x-dependence of the parameter c, which is perhaps a result of 
higher order corrections. A linear in x fit works very well, and the intercept represents 
a realistic estimate for c, 

c= 1.36(18) ■ (27r). (5.25) 

This is basically the same as the result of the x=0.01 point. 

The large value of the coefficient c explains why the heavy quark expansion does not 
work for the values of x considered in Fig. |^ along the dimensional reduction curve 
y = Vid^sdix). The adjoint "quark" can be considered heavy only if ^/y > c/{2tc), 
which gives x < 8/(9c^) ~ 0.01. 



21 



y = yjx) 32 



m/g. 



0.2 



0.1 



0.0 



-0.1 



2.0 



1.0 



0.0 



1 ' 

• 

♦ 

♦ 

♦ 

• 

• 


1 


• 

♦ 

1 


X 

c 

1 


1 ' 

X 

f ^ m(h) 

1 


1 

1 



0.00 0.10 0.20 

X 



0.30 



Figure 7: The order parameter (y4Qy4Q)/25f| and the scalar and vector masses m{AQ), 
m{h) in the symmetric phase along y = yAd~*3d{x) ■ When x < Xc ~ 0.235, the symmet- 
ric phase is metastable. The x < Xc -points have been calculated with (3g = 8, x > Xc 
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6 The physical phase of the 3d theory 

The previous simulation results were entirely for the 3d effective SU(2)+adjoint Higgs 
theory as such. We can now return to the main question: what is their relevance for 
the 4d finite temperature N = 2 QCD? 

According to the 2-loop dimensional reduction in Sec. 2, the physical 4d theories, 
defined by the values of T/A^, N = 2, Nj, lie on the straight lines yid^3d{x) (eq.( |2.27| )) 
plotted in Fig. ^. Consider first the quarkless case, Nf = 0. One observes that for 
almost all the range x < Xend ~ 0.30, one has y4:d->3d{x) < ydx) as already suggested by 
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Figure 8: Scalar, vector and glueball bound state masses for x — 0.01, 0.04, 0.10. 
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the 2-loop computations in Fig. |l|. Thus superficially, the finite temperature 4d theory 
corresponds to the broken phase of the 3d theory. This was the conclusion reached in 
(the status of which is reviewed in |[40|| ). 



However, this conclusion cannot be right. The 3d broken phase is a remnant of the 4d 
Z(2)-symmetry as discussed in Sec. 4, and the derivation of the effective theory is not 
reliable there. Instead, the physical phase must be the (in the 4d language, equivalent) 
symmetric phase where the 3d theory can be trusted. Indeed, at y = y4d->3d{x) the 
symmetric phase of the 3d theory is metastable in the whole region of x, and all 
the observables can be measured within the metastable phase, see Fig. ^ That the 
symmetric phase is the physical phase of the 3d theory was the conclusion reached 
in [Q, as well; however, there the 3d symmetric phase was seen to be stable. The 
difference is due to increased accuracy in relating continuum and lattice results. If one 
would be able to construct a 3d theory which fully respects the Z(A^) symmetry, then 
the metastable phase should be promoted to a stable phase in coexistence with its Z(2) 
companion. 

What happens when one goes to smaller T (larger x), towards T^? Then the 3d theory 
is becoming less reliable even in the symmetric phase, for several reasons. First, the 
perturbative expansion for the parameters of the 3d theory is becoming less reliable, 
see eq. ( |2.29|) . Moreover, the higher order operators are becoming more important 
since the mass hierarchy between the scales 5^3, vrT is lost due to increasing g'^{T) at 
T/A^ < a few, see Sec. 2. Indeed, the masses of the excitations in the 3d theory are 
~ g^, see Fig. |^, so they need then no longer be the dominant infrared excitations. 

In Fig. the order parameter and the masses have been plotted in the stable phase 
corresponding to y = yu^zd{,x). One sees that at a; ~ 0.235 the symmetric phase 
becomes stable. But this corresponds to T/Aj^ ~ 0.9, where the 3d theory is no longer 
reliable and the QCD phase transition has already taken place. Thus this "transition" 
is not related to the 4d QCD transition. 

Consider then the endpoint value, Xend ~ 0.30 (see Fig. H). This corresponds to 
T/Aj^ ^ 0.6 and is again below the value Td = 1.23(11) determined ||32[ for 
the SU(2) phase transition with 4d simulations. Yet the difference is not that large, 
and it is tempting to conjecture that if one had a 3d theory fully respecting the Z(A^) 
symmetry, then the corresponding endpoint might represent the QCD phase transition. 

Consider finally the Nf dependence. Fermions do not appear as dynamical fields 
in the effective theory, and their only effect is through the parametric dependence on 
Nf. At very small x (very large T) the curves in Fig. ^ imply that the system lies 
in the symmetric phase near Aq = 0. This result is in agreement with the fact that 
quarks break the Z(2) symmetry: the second minimum becomes metastable for Nf = 1 
and disappears for Nf > 3. At very large Nf, on the other hand, the perturbative 
expansion for the curves y = y^^sdix) gets worse. 
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broken phase, for x > 0.235 the symmetric phase (see Fig. ^). As discussed in Sec. 6, 
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7 Discussion 

Dimensional reduction allows to construct an effective field theory for high temperature 
QCD, which is reliable as long as -C vrT, implying T> a few x T^. However, the 
effective super-renormalizable 3d gauge+adjoint Higgs theory does not describe the 
confinement-deconfinement phase transition. 

The usual "power counting" picture of correlation lengths in high temperature QCD 
says that the longest scale is related to the magnetic sector of the theory and is of 
order {g'^T)~^ . The shorter scale, ~ {gT)^^, is associated with Debye screening. Our 
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study shows that this picture is correct only at extremely high temperatures. First, 
the purely magnetic scale comes from the lowest glueball mass in a pure SU(2) theory, 
rriG ~ 25f|. This should be compared with the leading order Debye mass mf, ~ ^g'^T^. 
The requirement m^) > tells that only at x < 0.005 are the "magnetic" effects 
numerically smaller than the "electric" ones. A similar estimate follows from the 
requirement that the non-perturbative corrections to the Debye screening mass are 
on the level of, say, 50 % of the tree-level value. In terms of temperature, this gives 
ridiculously large numbers. For example, for pure SU(2) this means T > lO^^Aj^jg. 

Thus, in the physically most interesting region somewhere above the critical temper- 
ature, the "naive" picture is wrong. The longest correlation length corresponds to the 
0^^ bound state Aq — Aq (in 3d language). In terms of 4d variables, this correlation 
length can be found from the analysis of Polyakov lines, 

(L(x)Lt(0))~exp[-m(Ag)|x|]/|x|, (7.1) 

where 

L(x)=Trn(x), =Vexp[ig [ drAoir,^)]. (7.2) 

Jo 

The second largest correlation length is associated with Debye screening. In 3d 
language, it is related to the operator hi. The corresponding 4d operators may be 
found in |37|. The non-perturbative Debye screening mass is about a factor 3 larger 
than the lowest order estimate up to temperatures T ~ 200Aj^. The fact that the 
Debye mass is non-perturbative does not allow a reliable integration out of the Aq field 
for the construction of the simplest effective field theory, containing only the scale g^, 
unless the temperature of the system is extremely large. 

Finally, the pure static glue sector has an even larger mass scale. In pure SU(2) and 
SU(2)-|-fundamental Higgs theories, rriG — ^-Qg^ in the present theory we have 
measured mc to be almost the same. 

It remains to be seen whether this modification of the standard picture of the high- 
temperature gauge theories has applications in the cosmological discussion of the quark- 
hadron phase transition or in heavy ion collisions. 
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Appendix A 

Some useful group theoretical relations for the SU(N) generators T°j in the fundamental 
representation are: 

1 1 i 

TaTf) 2,N ~^ '^dabsTs ~\~ "^fabsTgi 

TrTaTf, = -5a6, 
T^T^TaT^Tc = -{dabc + 'ifabc), 
T^TaTbTcTd = —^SabScd + -[dabsdcds — fabsfcds + Kdabs feds + fabsdcds)]- (A.l) 

4iV o 

Defining the adjoint representation by F^^ = —ifabe, one correspondingly gets 

Tr -F^Q-F^b fpaqfqbp -^^abi 

N 

FaFbFe ifpaqfqbrfrep i~^fabci 

''^^ FaFbFcFd fpaqfqbrfrcsfsdp 

N 

= SabScd + SadSbc + -^{dabsdcds — dacsdbds + dadsdbcs) 

1 

= -^{'^^ab^cd + ^ac^bd + ^ad^bc) + —{dabsdcds — fabsfcds)- (A. 2) 



Eq. ( A.l ) implies that for an adjoint vector Aa,A = AaTa, 



Id.abAaAbdscdAcAd = Tr - ^TtAY, (A.3) 
o iV 

SO that the matrix combination on the RHS isolates the 4-point coupling proportional 
to dsabdscd- For both SU(2) and SU(3), 

TrA^ = ^(TrA')^ (A.4) 
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